Unification of inflation, dark energy, and dark matter within the Salam-Sezgin 

cosmological model 



Alfredo B. Henrique^ 

Centra MulUdisciplinar de Astrofisica - CENTRA and Departamento de Fisica, 
Instituto Superior Tecnico, UTL, Av. Rovisco Pais, 104.9-001 Lisboa, Portugal 

Robertus Potting 
Centro Multidisciplinar de Astrofisica - CENTRA and Departamento de Fisica, 
FCT, Universidade do Algarve, Campus de Cambelas, 8005-139 Faro, Portugal 

Paulo M. Sii 

Departamento de Fisica, FCT, Universidade do Algarve, Campus de Cambelas, 8005-139 Faro, Portugal 

(Dated: March 11, 2009) 

We investigate a cosmological model, based on the Salam-Sezgin six-dimensional supergravity 
theory and on previous work by Anchordoqui, Goldberg, Nawata, and Nuiiez. Assuming a period 
of warm inflation, we show that it is possible to extend the evolution of the model back in time, to 
include the inflationary period, thus unifying inflation, dark matter, and dark energy within a single 
framework. Like the previous authors, we were not able to obtain the full dark matter content of 
the Universe from the Salam-Sezgin scalar fields. However, even if only partially successful, this 
work shows that present-day theories, based on superstrings and supergravity, may eventually lead 
to a comprehensive modelling of the evolution of the Universe. We find that the gravitational-wave 
spectrum of the model has a non-constant negative slope in the frequency range (10~^^ — 10^) rad/s, 
and that, unlike standard (cold) inflation models, it shows no structure in the MHz/GHz range of 
frequencies. 



I. INTRODUCTION 

We still are a long way from having a theoretically 
sound unified model explaining both the inflationary 
epoch and the standard era of expansion. The problem 
has become even more complex since the discovery of the 
dark energy content and the resulting present accelerated 
expansion of the Universe P, [I]- Recently, Anchordoqui, 
Goldberg, Nawata, and Nufiez proposed a very inter- 
esting model based on the Salam-Sezgin six-dimensional 
supergravity model [3| , whose spontaneous compactifica- 
tion from six to four dimensions gives rise to a potential, 
henceforth called the Salam-Sezgin potential. When this 
model is lifted to M-theory, the internal space is found 
to be non-compact Q , circumventing existing no-go theo- 
rems and allowing the potential to be positive. The work 
of Anchordoqui et al. Q is limited to the investigation 
of the epochs subsequent to primordial nucleosynthesis. 
The authors derived a solution of the Einstein field equa- 
tions in qualitative agreement with the observations of 
an accelerating universe, where the dark energy comes 
from a scalar field slowly rolling down the (exponential) 
Salam-Sezgin potential. Another important result was 
the finding that the model also contains a cold dark mat- 
ter component {jp = 0), coming from another scalar field 
interacting with the dark energy field, whose particles, as 



'Electronic address: 'alfredo@fisica.is t.utl.pt | 
^Electronic address: rpottingOualg.pt] 
J Electronic address: |pmsa@ualg.pt| 



a result of this interaction, have an effective mass slowly 
evolving with the dark-energy scalar field. Unfortunately, 
dark matter from the Salam-Sezgin scalar field only ac- 
counts for up to about 6% of the total matter content of 
the universe. 

In the present paper we extend this model back in time. 
It has been known for a long time 0] that an exponential 
potential of the generic form V = Vq exp(— v87rG7(/)) is 
able to support infiation only if 7 < The Salam- 

Sezgin potential fails to obey this requirement, having 
7 = ^/2. But, if we assume a period of warm infia- 
tion 0; it is possible to show that the Salam-Sezgin 
potential is able to drive enough inflation. We obtain, 
therefore, within the Salam-Sezgin cosmological model, 
a triple unification of inflation, dark matter, and dark 
energy. As in Ref. d, we were not able to obtain the 
full dark matter content of the Universe from the Salam- 
Sezgin scalar fields, the rest of the dark matter and the 
baryonic matter having to be put in by hand. Even if only 
partially successful, we believe that our analysis shows 
that present day theories, based on supergravity and su- 
perstrings, may eventually be successful in providing a 
comprehensive description of the evolution of the Uni- 
verse. 

During the period of warm infiation, energy is contin- 
uously transferred from the Salam-Sezgin scalar fields to 
a radiation fluid due to the presence of dissipative terms 
in the energy conservation equations (for a review see 
Q and references therein). For illustrative purposes we 
will assume a simple phenomenological dissipative term, 
following the suggestion of Yokoyama and Maeda Q . In 
this work we will not attempt to justify the occurrence 
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of warm inflation, which should be due to the coupHnffof 
the scalar fields to a heat bath of background fields [l^l ■ 
We just mention that plenty of such fields are present if 
the model is assumed to be derived from string theory. 

In the next section we briefly review the Salam-Sezgin 
model. In Sect. IIIII we derive the equations defining the 
model and introduce the details of how, following the in- 
flationary period, we calculate the amount of dark energy 
and dark matter. This is followed by a section on the nu- 
merical simulations of the evolution of the universe, since 
the beginning of the inflationary period till the present 
time, and the analysis of the results obtained. Section fVl 
is dedicated to the calculation of the full energy spectrum 
of the gravitational waves generated during the expan- 
sion of the universe, as it would be measured today by 
an ideal detector. Finally, we end with the summary and 
the conclusions. 



where {t, x) denotes a local coordinate system on M4 and 
Tc is the compactification radius of 82- The scalar fleld 
(j) is taken to depend only on the point of A/4, i.e., 4> = 
(j){t, x). The gauge field Am is excited on S2 and is of the 
form 



6cosi?, 



(4) 



which satisfies the Maxwell equations (obtained by vary- 
ing Am) and yields the field strength 



p2 



(5) 



Finally Bmp is taken to vanish. It follows that also 
Gmnp = 0. 

The Ricci scalar can now be written as \\!M. 



R[M] = R[Mi] + e-^fR[S2] - 40/ - 6(9^/)2 



(6) 



II. THE SALAM-SEZGIN POTENTIAL 



where i?[Af], i?[Af4], and R[S2] denote the Ricci scalars 
of the manifolds M , Mi, and S2, respectively. It follows 



We begin with a review of the Salam-Sezgin model, 
following the notation of Anchordoqui et al. Q. The 
Salam-Sezgin model represents a six-dimensional super- 
gravity[3l- It has been shown by Cvetic, Gibbons, and 
Pope that this model can be uplifted consistently to 
type I superstring theory. Its bosonic part is given by 
the action 



R - K (clMa) -e 



-K e Pmn '^MNP 



(1) 



Here, = detgjvflv, R is the Ricci scalar of gMN, c 
is a scalar field, Fmn = 9[M^Ar], Gmnp = OimBnp] + 



where An \& a gauge field and Bj^p is the 



i^^[mFnp\ , 

Kalb-Ramond field; capital Latin indices run from to 
5. Introducing Gg = 2k^ and ^ = 4g^ and rescaling 
(h = —Ka leads to 



S 



1 



(fx 



MN 



R- 

Gq 



idMcf'? 



^MNP 



(2) 



Here the length dimensions are: [Gg] = L^, [£,] = L^, 



L- 



= [gU] = 1, [Al,] = L-\ and [Fj,^] = [G^ 



MNP\ 



Salam and Sezgin [J showed that Lagrangian (P) can 
be compactified on the direct product M = Af4 x 5*2, 
where A'/4 = (Minkowski)4. It was shown by Gibbons, 
Giiven, and Pope [ll[ that this is the unique ground state 
among all non-singular solutions with a four-dimensional 
Poincare, de Sitter or anti-de Sitter symmetry. The 
monopole configuration they found has the metric on 
locally of the form 



R[S2] = 2/rl 



(7) 



and = e^^r'^^/gl, where 34 = det g^i, (with Greek 
indices running from to 3). We now define the gravita- 
tional constant in four dimensions as 



G4 = — — — 



Gfi 



327r2r2 ' 



(8) 



It follows that by using the Salam-Sezgin monopole con- 
figuration we can re-write the action in Eq. ^ as follows 



S = 



1 



167rG4 



x^gi< e 



J_p2/H 



R[Mi 



Gry 



-2/ 



Gfi 



e ■' 



2(9^/)' 
(9) 



Using the rescaling of the metric g^^, = e ' g^i, and ^/gl = 
e^-^ y/g4, the model is taken into the Einstein conformal 
frame where the action Eq. ^ takes the form 



S = 



1 



167rG4 



a;V.94 



Gfi 



-2f+4> , 



Gry 



-6/- 



-4/ 



(10) 



The four-dimensional Lagrangian is then 



L 



WttG 



[i? -4(9^/)2- (9^0)2 _y(/,^)], (11) 



with 



-2/+0 



Gy_ 



-If 



Gfi 



(12) 



dsl = dsi{t, xf -I- e'^^'-*'^^rl{d'd^ + sin^ Mip'^), (3) where we have written g = g4,R = R[g4] and G = G4 
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Following Anchordoqui et al. Q, we now perform the 
field redefinition 



x={(t> + 2/)/\/f67rG, y = (0 - 2/) /\/l67rG. (13) 
The Lagrangian becomes 



R 



].{dxf-].[dyf-V{x,y) 



167rG 

with the Salam-Sezgin potential given by 

V{x, y) ^ Aie'^y (l - 2yl2e-"^ + Age-^"^ 
where we have introduced the notation 



5127r3r2 



327r2 



4^2 



^3 = 



1024:TT% 



(14) 



(15) 



(16) 



and a = VTGttG — Ay/Tr/nip. 



III. THE EVOLUTION MODEL 

In our analysis of the Salam-Sezgin cosmological 
model, we will divide the evolution of the universe in 
two stages. The first stage starts at the beginning of the 
inflationary period and extends well into the radiation- 
dominated epoch. During this first stage of evolution, 
energy is continuously transferred from the scalar fields 
X and y to a radiation fluid due to the presence of dis- 
sipative terms in the energy conservation equations. In- 
flation, which we assume to be of the warm type 0, is 
driven by the field x. When the energy density of ra- 
diation is of the order (lO^^'GeV)^ (see Sect. |lV]for de- 
tails), we assume that the energy transfer from the fields 
X and y to radiation ceases. This marks the end of the 
first stage of evolution. During the second stage, which 
extends up to the present epoch, the z-field oscillates 
around the minimum of the potential, behaving as cold 
dark matter with varying mass, as explained below. The 
y-field, which behaves like dark energy, is practically con- 
stant during most of the second stage of evolution and 
becomes dominant in a recent epoch, in agreement with 
observations. In what follows we present the relevant 
equations for both stages of evolution and show how an 
unification of inflation, dark matter, and dark energy can 
be achieved within the Salam-Sezgin cosmological model. 

For a flat universe, described by the Friedmann- 
Robertson- Walker metric. 



(17) 



where a{t) is the scale factor and dy^ is the metric of the 
three-dimensional Euclidean space, the Einstein equa- 
tions are given by 
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Y + y + Vix, y) + pr 



(18) 
(19) 



where V{x, y) is the Salam-Sezgin potential given by 
Eq. (jlSp . Pr is the energy density of a radiation fluid 
with equation of state pr ~ Pr/i, and a dot denotes a 
derivative with respect to the comoving time t. 

Let us assume that, during the flrst stage of evolution, 
energy is continuously transferred from the scalar flelds 
X and y to the radiation fluid pr due to the presence of 
dissipative terms. In this case, the energy conservation 
equations for the scalar fields and the radiation fluid are 



i^-x 
a 



&V 
dx 

.. oO. dV 

y + 3-y+ — 
a oy 



Pr 



A " 

A-pr 

a 



(20) 
(21) 
(22) 



where F^, and Fj, are the dissipative coefficients. 

The form of the dissipative coefficients has been dis- 
cussed in literature and can be, in realistic models, quite 
complicated. For illustrative purposes we will consider 
in this article simple phenomenological dissipative coef- 
ficients, first introduced by Yokoyama and Maeda @, of 
the form 



^x{x,y) = .^1 



^v{x,y) = fy 



l aw 

dx'^ ' 



dy"^ 



(23) 
(24) 



where fx and fy are free parameters. For an appropriate 
choice of these parameters, a significant amount of en- 
ergy is transferred from the scalar fields x and y to the 
radiation fiuid and, as a consequence, the energy density 
of radiation is not diluted during inflation, as in the stan- 
dard (cold) inflationary models, and its influence in the 
evolution of the universe is such that enough inflation oc- 
curs, despite the steepness of the Salam-Sezgin potential. 

To be general, let us assume that, immediately prior to 
the inflationary phase, the universe may have, in addition 
to the energy contained in the x and y scalar fields, an 
important contribution from the radiation fluid. Even- 
tually, the potential energy V{x,y) becomes dominant 
and warm inflation begins. Initially, the energy density 
of radiation, p^, decreases faster than the energy den- 
sity of the fields x and y, pxy = x^ /2 + jl -I- V{x,y^. 
However, because of the energy transfer due to the dis- 
sipative terms, pr quickly reaches a state during which 
it decreases slower than p^y At a certain point, p^ be- 
comes greater than p^y and a smooth transition from 
warm inflation to a radiation-dominated universe takes 
place [ist . 

During inflation the x-field slowly rolls down the po- 
tential y(x,y), approaching its minimum at Xmin = 
(1/a) ln(A3/^2)- During the inflationary period the y- 
field does not play a significant role in the evolution of 
the universe; it simply slowly rolls down the exponential 
(in the y direction) potential. 
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Note that the potential (jlSp can be expanded near its 
minimum, yielding 



V{x, y) = Ae^y + ]^MI{x - x^,,^f + (25) 



where 



^3 



(26) 



is the time-dependent mass of the scalar field x and 

(fe^e - 1) . (27) 



The first term on the right-hand-side of Eq. (|25|) behaves 
like dark energy and, at the present epoch, dominates the 
dynamics of evolution of the universe. In order to achieve 
agreement with observations, the constant A should be 
very small, of the order of lO^^^^mp (see Sect. |IV] for 
details). This extremely small, but different from zero, 
value means that supersymmetry should to be broken 
within the Salam-Sezgin cosmological model 3. 

After the end of the inflationary period, energy transfer 
from the x and y fields to the radiation fluid continues. 
As the universe expands, its temperature drops and even- 
tually the interaction between the scalar fields and the 
thermal bath becomes negligible. At this point, we turn 
the dissipative coefficients and Ty to zero. This marks 
the end of the first stage of evolution. It is necessary to 
emphasize here, that the residual amount of energy in the 
x-field at the end of the first stage of evolution should be 
low, in order to make sure that the universe undergoes 
a long enough radiation-dominated period, ending well 
after nucleosynthesis [3, [TB| . Our numerical simulations 
(see next section) show that the transition between the 
first and second stage of evolution occurs at temperatures 
of the order of 10^^ GeV and that the scalar field x is, at 
that time, of the order of — 10~^^mp. 

Due to the absence of the dissipative terms, the scalar 
field X oscillates rapidly around its minimum during the 
second stage of evolution. Since in the expanded poten- 
tial (|25p the dominant term is quadratic, the oscillating 
cc-field behaves like matter with equation of state p = 
[l6t . This pressureless matter behaves like cold dark mat- 
ter. Its energy density is proportional to the y-depcndent 
mass of the x-field (which is practically constant during 
most of the second stage of evolution) and decreases as 



\ a 



(28) 



where c is the speed of light, a — 4\/7rG/c, and ao = a (to) 
is the value of the scale factor today^. The constant C 



^ In the first stage of evolution we use the natural system of units, 
with mp = 1/VG = 1.22 X 10^^ GeV, while in the second stage 
we use the International Systems of Units. 



should be chosen in order to guarantee continuity of the 
solution at the transition from the first to the second 
stage of evolution, namely, 



C 



p2 p-aVT/2 



Pr,0 



3/4 



(29) 



where the subscript T refers to quantities evaluated at 
the end of the first stage of evolution and we have taken 
into account that 



Pr 



Pr.Q — 

v a 



(30) 



where p^.o = Pr{to) ~ 4.6 x 10^'^^ kg/m"^ is the density 
of radiation observed today. 

Within the Salam-Sezgin cosmological model, the x- 
field accounts for just a fraction of the total matter con- 
tent of the universe 0]. Therefore, in order to be in 
agreement with observational data, we introduce an ad- 
ditional pressureless matter component, which accounts 
for other types of dark matter and for the usual (bary- 
onic) matter. For this component, the energy conserva- 
tion equation yields 



B fao 
Pni = — 1 — 



(31) 



where J3 is a constant. 

The evolution equations for the y-field and the scale 
factor a are now 



y = 

a 
a 



a 2 V a / 

2f - 2Ae°'y +(b + Ce"y^'^ 



(32) 



a 
12 



+ 2p,oc^ (^) 
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+ Pr 



Ae^y + + Ce"^/2^ 



a J 



2 /aoV 

It) 



(33) 



(34) 



In the above equations, the constants A and B should 
be chosen consistently with observations, namely. 



Ae^y" 



2 



PoEfiC 



(35) 
(36) 



where /9de,o = 6.9 x 10 kg/m"^ and Pm,o = 2.6 x 
10^27 kg/ni^. To the values of 

Pde,Oj Pm,Oj and pr^o used 
in this article, correspond a value of the Hubble constant 
Ho = 71 kms-^Mpc~^ 

Following Ref. [3], during the second stage of evolution, 
we use, instead of the comoving time t, a new variable u, 
defined as 



■ln(l + z). 



(37) 
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where z = oo/a — 1 is the redshift. 

Using this new variable, the energy conservation equa- 
tion for the scalar field y can be rewritten as 



and 



a „ / a 
- + 2 - 
a \ a 



aAe° 




(38) 



where the subscript u denotes a derivative with respect 
to u; d/a and {a/a)'^ are functions of u, y and given 
by 



a a 
a ~ 12 



= 2a' 



|4a2 Ae^y + {b + Ce"^/^^ g-^" 

yl[a^yl - 12)-^ +2Ae'^y 
[b + Ce"«/2^ - 2p,,oc'e-*"} , (39) 

Ae°'y +{b + Ce"^/2^ e'^" 

{12 -a^yl)-\ (40) 



The set of equations Eqs. (fT8|) - ([22)) . for the first stage 
of evolution, and p8 |) -(|40 |) . for the second stage, will 
be solved numerically in the next section. The parame- 
ters fx and fy will be chosen such that enough inflation 
(65 or more e-folds) is achieved during the first stage 
of evolution. The time of transition between the two 
stages of evolution will be chosen such that the radiation- 
dominated epoch taking place after inflation ends well af- 
ter nucleosynthesis. The constant C will be chosen such 
that the solution is continuous at the transition between 
the first and second stages of evolution and the constants 
A and B will be chosen in order to guarantee that the en- 
ergy densities of radiation, matter (baryonic plus dark), 
and dark energy at the present time are in agreement 
with observations. 

In section IIVI our numerical results will be presented 
in terms of the density parameters for radiation, mat- 
ter, and dark energy, the equation-of-state parameter for 
dark energy, and the deceleration parameter defined as, 
respectively. 



Pr 
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B + Ce"y/^] ( ^ 



^ + Ae^y 
2 \ a 



(41) 
(42) 
(43) 



where p^ = SH^/iSnG) = &{a / af / {a^ c^) is the critical 
density, 



w 



^ _Py_ 



yl{^)~2Ae-y 



2Ae°'y 



(44) 



aa 
72- 



(45) 



To finish this section, let us point out that a triple 
unification of infiation, dark matter, and dark energy 
was proposed recently by Liddle, Pahud and Urena- 
Lopez [15|, using a single massive scalar field (f> with 
potential V{(j)) = Vo + ^m^^^. Note that this poten- 
tial is similar to the expanded Salam-Sczgin potential 
given by Eq. ((25)) . the main difference being the fact that 
the term behaving like dark energy is constant and not 
time-dependent as in our case. Within the unification 
scenario proposed in Ref. pj| . a residual inflaton field 
(j) survives preheating and oscillates rapidly around the 
minimum of the potential, behaving as cold dark mat- 
ter. However, similarly to our case, the amplitude of 
the scalar field oscillations are too high. They should be 
drastically reduced in order to allow for a long enough 
radiation-dominated epoch encompassing primordial nu- 
cleosynthesis. Exploiting the uncertainty in the cosmo- 
logical evolution between the end of the inflationary pe- 
riod and nucleosynthesis, the authors of Ref. [3l consid- 
ered a brief period of thermal infiation following preheat- 
ing, in order to reduce the amplitude of the scalar field 
oscillations to the desired level. In this work, we man- 
aged to reduce of the amplitude of the x-field oscillations 
after infiation within the context of the Salam-Sezgin cos- 
mological model. As already explained above, this was 
accomplished by extending the energy transfer from the 
Salam-Sezgin scalar fields to the radiation fluid until well 
into the radiation-dominated epoch. 



IV. NUMERICAL SIMULATIONS 

Let us now solve numerically the set of equations (fT8|) - 
and (|38p - (PDl) . As initial conditions, at the be- 
ginning of the inflationary period, we choose x{ti) = 
— O.lmp and y{ti) = 0.1 mp. Note that enough infla- 
tion (65 e-folds or more) can be achieved for such initial 
amplitudes of the inflaton field x, below the Planck mass, 
provided that fx is chosen adequately. In what follows 
we choose the phenomenological dissipative parameters 
to be fx = fy = 275. In order to achieve agreement 
with observations, the constant A, introduced in Eq. (|27p . 



should be of the order of 10 ^^^rn^. We choose, there- 
fore, Ai = 10-12 A2^l- 10-"", and ^3 = 1. This 
choice of Ai also fixes the energy scale of infiation to be 
E^ijjf = V{xi,yrf/'^ 10^'' GeV. We assume that, at the 
beginning of the inflationary period, the kinetic terms of 
the scalar fields, as well as the energy density of radiation, 
are of the same order of magnitude of the potential V . 
We choose, therefore, x{ti) = y{ti) — 1.47 x lO^^mp and 
Pr{ti) = 4.33 X lQ~^'^m'^. Finally, we choose a{ti) = 1; 
a{ti) is fixed by the Eq. ^ to be 8.52 x lO-^^mp. 

For such values of the initial conditions and parame- 
ters, the infiationary period lasts for about 1.1 x 10^ tp 
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0.0 0.5 1.0 1.5 2.0 



FIG. 1: Time evolution of the energy densities of radiation pr 
(dashed curve) and of the scalar fields pxy = i^/2 + y^/2 + 
V{x,y) during the first stage of evolution. The smooth tran- 
sition from inflation to a radiation-dominated universe occurs 
at about l.f x 10**tp. 




0.0 0.5 1.0 1.5 2.0 

t/(10'*(p) 



FIG. 2: Time evolution of the scalar fields x (dashed curve) 
and y during the first stage of evolution. After the inflation- 
ary period y decreases very slowly, while the |a;| decreases 
exponentially. 



and the scale factor grows by a factor of about 10^*. Dur- 
ing inflation the condition > H holds, as expected 
in warm inflation. Furthermore, the dissipative coeffi- 
cient F-c is much greater than the Hubble parameter H , 
implying a strong dissipative warm- inflation regime Q. 
At the very beginning of the inflationary period, the en- 
ergy density of radiation decreases faster than the energy 
density of the scalar fields. However, because of the en- 
ergy transfer due to the dissipative terms, the decrease 
of the energy density of radiation slows down and, after 
a while, a smooth transition from warm inflation to a 
radiation-dominated universe takes place (see Fig.[T|). 

During inflation the x and y-fields roll down the Salam- 
Sezgin potential. At the end of the inflationary period 
their values are x « —9.30 x 10~^TOp and y « 2.16 x 
10~^mp. During the subsequent evolution, till the end 
of the first stage of evolution, y decreases very slowly, 
while l^l decreases exponentially (sec Fig. [5]). 
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u 

FIG. 3: Time evolution of the density parameters Qr (thin 
curve), fla (dashed curve), and fly (thick curve) during the 
second stage of evolution. In this case, which corresponds to a 
short duration of the first stage of evolution, 8.00 x 10* tp, the 
transition from a radiation to a matter-dominated universe 
occurs earlier in the history of the universe (but well after 
nucleosynthesis) and the density parameter for dark energy 
becomes non-negligible already at a redshift z ~ 10®. The 
x-field dark matter amounts to about 10% of the total matter 
content of the universe at the present time uo = 0. 



As already mentioned above, the residual x-field oscil- 
lates rapidly around the minimum of the potential during 
the second stage of evolution in which fx ~ fy ~ 0, be- 
having like cold dark matter. In order to achieve a long 
radiation-dominated period of evolution, the initial am- 
plitude of the a;-field oscillations should be much smaller 
than the Planck mass Our numerical simulations 

show that a long enough radiation-dominated epoch re- 
quires, at the time tr of transition between the first and 
second stages of evolution, that Xt ~ — 10~^^TOp. There- 
fore, the energy transfer from the scalar fields x and y to 
the radiation fluid should continue after the end of the 
inflationary period, allowing \x\ to decrease from about 
10"'^ nip, at the end of inflation, to about 10~^^mp, at the 
time of transition between the first and second stages of 
evolution. 

The value of should be chosen very carefully. If 
this value is too small, the amplitude of the scalar field x 
does not decrease enough before it starts to oscillate and, 
consequently, the transition from a radiation to a matter- 
dominated universe takes place too soon. Another con- 
sequence of a shorter first stage of evolution is a rise of 
the energy density of the y-field (dark energy) at redshift 
z ^ 1. In Fig. [21 the evolution of radiation, matter (bary- 
onic plus dark) , and dark energy during the second stage 
of evolution is shown for = 8.00 x 10^ tp. In this case, 
Xj = —1.58 x 10~^^TOp and p^^o = 9.7 x IO^^Pm.o- On 
the other hand, if the value of is too big, the amplitude 
of the x-field decreases too much and its contribution to 
the total matter content, at the present time, becomes 
neghgible. For instance, in the case tr = 9.50 x 10^ tp we 
have Xt = —4.20 x 10~^^ mp and px,o = 0.2 x 10~^ Pm,o- 

Let us now analyze in more detail the solution obtained 
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FIG. 4: Time evolution of the j/-field dark energy. This field, 
which is practically constant during most of the second stage 
of evolution, became dominant only in a recent past. 
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FIG. 5: Time evolution of the equation-of-state parameter for 
dark energy, Wy, during the second stage of evolution. At the 
present time, ito — 0, the value of this parameter is about 
—0.62. In the future, it will approach the value —1/3. 



for = 8.90 X 10^ tp. In this case, at the time of transi- 
tion between the first and second stages of evolution, we 
haveoT = 9.85 xlO^^ Xt = -4.45xl0-i'^mp, = 2.09 x 
IQ-^nip, p^^r = 1.16xlO"''3^4 ^ 5.96 X 10^3 ]Jg/m^ and 

Pr^T = 4.89 X 10-20 = 2.52 x 10^"^ kg/in . Taking into 
account that pr,o = 8.92 x lO'^^^ m;? = 4.6x lO'^^ kg/m^, 
we obtain Ut = (1/4) ln(pr.o/Pr,T) = —62.02, and, from 
Eq. C = 1.68x10-124^4 ^ 7.80xl0-"kgm-is~2. 

The constants A and B are chosen in order to sat- 
isfy the observational constraints ([55]) and ([M]) . namely, 
A = 3.16 X 10-122 = 1.47 X 10-*^kgm"is-2 and 
B = 4.73 X 10-124 mf = 2.19 x lO-i" kg.m"is-2. In this 
case, we obtain that the x-field dark matter contributes 
about 6% to the total matter content of the universe at 
the present time, p^^o = 6.15 x 10-2 p^ g. The time evo- 
lution of the quantities y, Wy, q, O^, and fly, defined 
in Eqs. (|4T|) - (|45)l . is shown in Figs. [4H2l for the case under 
consideration. 

To finish this section, let us comment on the evolu- 
tion of the scale factor today and in the future. As can 



FIG. 6: Time evolution of the deceleration parameter, q, dur- 
ing the second stage of evolution. This parameter has its 
minimum (negative) value around the present time and ap- 
proaches zero in the future. 
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FIG. 7: Time evolution of the density parameters Qr (thin 
curve), (dashed curve), and fly (thick curve) during the 
second stage of evolution. In this case, which corresponds 
to a duration of the first stage of evolution of 8.90 x 10* fp, 
the density parameter for dark energy becomes non-negligible 
only recently. At the present time. Mo = 0, the density param- 
eters for radiation, flr.o ~ 4.8 x 10-^, matter, flufi = 0.27, 
and dark energy, fly.o ~ 0.73, are in good agreement with 
observational data. The x-field dark matter contributes only 
about 6% to the total matter content of the universe. 

be seen in Fig. O the equation-of-state parameter for 
dark matter, Wy, is today about —0.62, implying that 
a{t) oc i^'^^, if the matter and radiation contributions 
are neglected. However, there is today a non-negligible 
contribution from matter {flufi = 0.27), which if taken 
alone would imply a{t) oc ^2/3. Taking into account both 
contributions, from dark energy and from matter, wc ob- 
tain that today the expansion of the universe is slightly 
accelerated. Naively it would seem that this acceleration 
should increase in the future, since the density parameter 
for matter will approach zero and dark energy will domi- 
nate the dynamics of evolution of the universe. However, 
at the same time as f^m — > 0, the cquation-of-statc param- 
eter for dark energy, Wy, approaches —1/3 (sec Fig. [5]). 
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FIG. 8: Time evolution of tlie scale factor a (thick curve). To 
guide the eye, a linear growth is also represented (dashed 
straight line). The transition to an accelerated expansion 
takes place at about (2 — 3) x 10^^ s. At the present time, 
t ~ A X 10^^ s, the expansion of the universe is slightly accel- 
erated. In the future, the scale factor growth approaches a 
linear behavior, a consequence of the fact that the equation- 
of-state parameter for dark energy approaches —1/3. 



that gives the number of gravitons. All quantities in 
the above expressions are evaluated at the present time, 
uq = 0. The squared Bogoliubov coefficient is given 
by = (X — y)^/4, where X and Y are continuous 
functions of time, determined during the first stage of 
evolution by the set of differential equations 



X 
Y 



I 



\ a / 
a 



2 /flo 



\ ^ a 
J a 



■ \ 2 

a 



(46) 



X, (47) 



where d/a and (a/a)^ are given by Eqs. p8)) and p^ . 
respectively, and during the second stage of evolution by 
the set of differential equations 



Xu — 



'ao\ Y 
a / a/a 



I a 
loq ao 
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X 

a/a' 



(48) 
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As a result, the deceleration parameter, q, approaches 
zero (see Fig. ^ and the scale factor grows linearly in 
the future. In Fig. [8] we show the evolution of the scale 
factor into the future, obtained by direct numerical in- 
tegration of Eqs. ([32|l and p4|) . using the same initial 
conditions and values of the parameters A, B, and C as 
in Figs. HHll 

V. GRAVITATIONAL WAVES 

In this section, we calculate the full gravitational-wave 
spectrum of the Salam-Sezgin cosmological model using 
the formalism of the continuous Bogoliubov cocfhcicnts. 
This formalism was first applied to particle production in 
an expanding universe by Parker |17| and later extended 
to the case of gravitons by Henriques [ll] and Moor- 
house, Henriques, and Mendes [l9| . It avoids, in a natural 
way, overproduction of gravitons of high frequencies [20l | . 
In recent years, the formalism of continuously evolving 
Bogoliubov coefficients was applied to the calculation of 
gravitational-wave spectra in several cosmological mod- 
els, revealing interesting features in the MHz/GHz range 
of frequencies due to the transition between the inflation- 
ary and the radiation-dominated eras 

Let us define the gravitational-wave spectral energy 
density parameter as 

^ ^ 1 dpcw ^ 8hG _^4^2 
Pcdlnuj Sttc^H'^ ' 

where Pq^^c^ and to are the energy density and angular 
frequency of the gravitational waves, respectively, pcC^ 
is the critical energy density of the universe, H is the 
Hubble parameter, and /3 is a Bogoliubov cocfhcient, such 



where a/ a and {a/a)^ are given by Eqs. ([55)1 and (|40p . 
respectively. The scale factor at the present time is given 
by ao = 0x6""^, where the subscript T stands for the 
time of transition between the first and second stages of 
evolution. 

The above differential equations are integrated numer- 
ically, from the beginning of the inflationary period till 
the present time. We assume that prior to the infla- 
tionary period the universe was radiation-dominated and, 
therefore, no gravitons were created, i.e., a'^[t) = 1 and 
0^{t) = for t < ti, where a is a Bogoliubov coefficient 
related to X and Y through ^ {X + Y f/i. The ap- 
propriate initial conditions for the functions X and Y are 
then X{ti) = 1 and Y{ti) = 1. The gravitational- wave 
angular frequency at present time, loq, is taken to vary 
from about 10~^^rad/s (corresponding to a wavelength 
equal, today, to the Hubble distance) to about 10®rad/s 
(corresponding to a wavelength equal to the Hubble dis- 
tance at the end of the inflationary period). 

The full gravitational-wave spectrum for the Salam- 
Sezgin cosmological model is shown in Fig. [HI It was 
calculated for the particular case of the model with = 
8.90 X 10^ tp, discussed in detail at the end of Sect. HVl 

For angular frequencies corresponding to the present 
size of the Hubble horizon, an upper limit on the 
gravitational-wave spectral energy density parameter can 
be derived from measurements of the cosmic microwave 
background radiation, narnely, < 1-4 x 10~^" for 

i^hor ~ 1.4 X 10~^^rad/s [2J]. The gravitational-wave 
spectrum of the Salam-Sezgin cosmological model satis- 
fies this constraint easily. 

In the intermediate range of frequencies, from about 
10~^^rad/s to about 10®rad/s, the gravitational-wave 
spectrum differs substantially from the ones obtained, 
with the same formalism, in Refs. [U, [12, , in that it 
has a non-constant slope. This is due to the fact that, in 
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FIG. 9: Gravitational wave spectrum of the Salam-Sezgin 
cosmological model. 



FIG. 10: Time evolution of the equation-of-state parameter, 
w = p/p, during the first stage of evolution. 



previous investigations, numerical evaluation of the Bo- 
goliubov coefficients started at the end of inflation and, 
moreover, it was assumed, for the purpose of determin- 
ing the initial conditions for the functions X and Y , that 
the inflationary period was either exponential or of the 
power-law type. This assumption was needed since exact 
analytical solutions of Eqs. (pS)) and (|T7)) are known only 
for these types of inflation. As a result, the gravitational- 
wave spectrum had, in the intermediate range of frequen- 
cies, zero slope [111 [2^ or constant negative slope psj . 
respectively. In this article, we have improved our cal- 
culation by numerically evaluating the Bogoliubov coeffi- 
cients from the very beginning of the inflationary period. 
This approach allow us to consider any type of inflation, 
not just exponential or power-law. This is rather con- 
venient for the cosmological model under consideration, 
since, during inflation, the equation-of-state parameter 
w = p/p changes smoothly from —1 to —1/3 (see Fig.irU]). 
i.e., inflation is neither exponential nor power-law. 

The shape of the spectrum in the intermediate range 
of frequencies can be understood as follows. It is known 
that, for power-law inflationary models, the spectrum has 
constant negative slope in this range of frequencies [25 1 
(see also Ref. [g^] for the derivation of this result using 
the formalism of continuous Bogoliubov coefficients). In 
the limit where power-law inflation approaches exponen- 
tial inflation (s +oo, where a{t) oc t^), the slope of 
the spectrum approaches zero. Now, let us regard the 
evolution of our cosmological model, during inflation, as 
a succession of short periods of time, with duration At, 
such that in each of them inflation can be considered to 
be of the power-law type with s — |(w4-l)~^ w constant. 
Gravitational waves that cross the Hubble horizon dur- 
ing the first of these periods, have today low frequencies. 
Since during this period the equation-of-state parameter 
is w = p/p « — 1 (see Fig. [TU]), corresponding to s » 1, 
the spectrum, for such low frequencies, is almost flat. 
During the second period of time, w is slightly higher 
and the spectrum has a small negative slope in the corre- 
sponding frequency range. Continuing with this process. 



we obtain a spectrum that is a succession of line segments 
with increasing (in modulus) slope as the frequency of the 
gravitational waves increases. Taking the limit At — *■ 0, 
we obtain, in the intermediate region of the spectrum, 
the shape shown in Fig. [51 

The gravitational-wave spectrum of the Salam-Sezgin 
cosmological model shows no structure in the MHz/GHz 
frequency range. Such a structure is expected in cosmo- 
logical models in which standard (cold) inflation is fol- 
lowed by a preheating and/or reheating stage f2ll. [2^. [23j . 
In such models, after the inflationary period the infla- 
ton field oscillates about the minimum of the its poten- 
tial, leading to particle production and an increase of 
the energy density of radiation. The structure of the 
gravitational- wave spectrum in the MHZ/GHz range of 
frequencies depends on the form of the infiationary po- 
tential near its minimum. More specifically, if the po- 
tential near its minimum is proportional to 0", then 
ilcw increases with the increase of n (23j . Contrarily to 
these standard (cold) inflationary cosmological models, 
we have assumed that, within the Salam-Sezgin cosmo- 
logical model, inflation is of the warm type. That is, 
the energy density of the x and y-fields is continuously 
transferred to the radiation fluid and, consequently, the 
energy density of radiation is not diluted during inflation. 
Because the energy density of radiation decreases slower 
than the energy density of the x and y-fields (see Fig. [T|) , 
the transition from inflation to the radiation-dominated 
era occurs smoothly, well before the x-field starts to os- 
cillate about the minimum of its potential. Therefore, 
the coherent oscillations of the z-ficld leave no imprint 
in gravitational- wave spectrum in the MHz/GHz range 
of frequencies. 



VI. CONCLUSIONS 

In this work we have extended the analysis by Anchor- 
doqui, Goldberg, Nawata and Nunez based on the Salam- 
Sezgin supergravity theory, by including the infiationary 
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period in the evolution. The Salam-Sezgin potential con- 
tains two scalar fields, x and y. Assuming sufficient dissi- 
pative coupling, it follows that a period of warm inflation 
takes place, driven by the x-field, during which energy 
is being transferred continuously from both scalars to a 
radiation fluid. Well after the end of the inflationary pe- 
riod, with the universe already dominated by the energy 
density of the radiation fluid, the dissipative coefficients 
are switched off and the x-field begins a period of quick 
oscillations around the minimum of the potential, behav- 
ing like cold dark matter with a varying mass. The y- 
field, which behaves like dark energy, becomes then prac- 
tically constant until very near the present epoch, when 
it begins to dominate the energy content of the universe. 
We have shown, therefore, that within the Salam-Sezgin 
cosmological model is possible to obtain a unified descrip- 
tion of inflation, dark matter, and dark energy. For an 
appropriate choice of the parameters of the model, we 
were able to account for the observed values of radiation, 
matter, and dark-energy density parameters. However, 
only approximately 6% of the total matter content of the 
universe was obtained from the i-field, while, according 
to observational data, one would expect dark matter to 



amount to about 80% of the total matter content. Even 
if not entirely successful in this respect, it is not unrea- 
sonable to expect that more sophisticated models, based 
on supergravity and superstring theories, could be found 
that may generate enough dark matter and dark energy 
in a fully unified description. 

We have also calculated the full gravitational-wave 
spectrum of the Salam-Sezgin cosmological model. Using 
the formalism of the continuous Bogoliubov coefficients, 
we have shown that this spectrum has a non-constant 
negative slope in the frequency range (10~^^ — 10^) rad/s, 
and that, contrarily to cosmological models in which 
standard (cold) inflation is followed by a preheating 
and/or reheating stage, the gravitational- wave spectrum 
shows no structure in the MHz/GHz range of frequencies. 
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